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O ' Abstract 

^: 

CD ■ Trajectory-based approaches to quantum mechanics include the 

de Broghe-Bohm interpretation and Nelson's stochastic interpreta- 
tion. It is shown that the usual route to establishing the validity of 
such interpretations, via a decomposition of the Schrodinger equation 
into a continuity equation and a modified Hamilton-Jacobi equation, 

^ I fails for some quantum states. A very simple example is provided 

by a quantum particle in a box, described by a wavefunction initially 

Q^. uniform over the interior of the box. For this example there is no 

corresponding continuity or modified Hamilton-Jacobi equation, and 
the spacetime dependence of the wavefunction has a known fractal 
structure. Examples with finite average energies are also constructed. 



1 Introduction 

The formalism of standard quantum mechanics is very different from that of 
classical mechanics, as are the generic phenomena described by each theory. 
This gives rise to a number of well known interpretational issues when one 
tries to integrate classical and quantum aspects of the world. The Copen- 
hagen interpretation of quantum mechanics remains foremost for most physi- 
cists in resolving such issues, and adequately explains the empirical content 
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of the standard formalism. However, a number of alternative interpretations 
exist, and can be valuable in providing (i) reasonably coherent pictures for 
thinking about fundamental quantum phenomena such as interference and 
entanglement; (ii) means for marrying the microscopic with the macroscopic 
(in contrast to the enforced separation specified by the Copenhagen inter- 
pretation); and (iii) a variety of starting points for extending or modifying 
the standard quantum formalism. 

One class of alternative interpretations is distinguished by retaining the 
classical concept of spacetime trajectories for quantum particles. Thus, for 
example, in the de Broglie-Bohm interpretation particles follow trajectories 
in configuration space determined by a guiding wave, providing an underlying 
deterministic (but nonlocal) picture of quantum evolution [H I2] • A second 
example is Nelson's stochastic interpretation, in which particles follow non- 
differentiable trajectories in configuration space determined by a stochastic 
generalisation of Newton's second law Ol^. 

An important claim made by proponents of such trajectory-based in- 
terpretations is that they reproduce all predictions of standard quantum 
mechanics, at least for the case of nonrelativistic particles moving under 
velocity-independent potentials [H El El • This general reproducibility is 
a necessary feature of any complete interpretation of the quantum formal- 
ism. Indeed, the dBB and stochastic interpretations appear capable of going 
beyond the standard formalism, as they have even been applied to non- 
normalisable wavefunctions lying outside the Hilbert space of possible quan- 
tum states (eg, sections 4.10, 4.11 of Ref. [2j, and page 88 of Ref. j3]). 

However, the aim of this note is to show that such trajectory-based inter- 
pretations in fact do not apply to all quantum states, and hence are formally 
incomplete. Moreover, since some of the states in question have finite aver- 
age energies, it is suggested that such interpretations may also be physically 
incomplete. Similar difficulties arise for any interpretation based on conti- 
nuity and modified Hamilton- Jacobi equations, including the hydrodynamic 
and exact uncertainty approaches to quantum mechanics |21ISlini- 

The results are based on a subtle property of the Schrodinger equation in 
the position representation, discussed in section 2. In particular, for certain 
states of quantum particles, this equation cannot be partitioned into separate 
terms involving spatial and temporal derivatives of the wavefunction respec- 
tively. A simple example is provided by a wavefunction initially uniform over 
the interior of some region and vanishing elsewhere (eg, a plane wave incident 
on a slit, or a particle confined to a box with maximal position entropy). 



The existence of such states formally arises as a consequence of the un- 
boundedness of the corresponding Haniiltonian operator - this operator can- 
not in fact be directly applied to a large class of wavef unctions, even though, 
paradoxically, these wavef unctions and their evolution are themselves per- 
fectly well defined. Further, as shown in section 3, such wavefunctions can 
have finite average energies, and hence can, in principle, be physically pre- 
pared from finite resources. 

An important consequence of unboundedness is that the Schrodinger 
equation cannot always be decomposed into a continuity equation and a 
modified Hamilton- Jacobi equation. As discussed in section 4, the nonexis- 
tence of such a decomposition for certain states leads to an incompleteness 
of trajectory-based interpretations for these states. 

The lack of continuity and modified Hamilton- Jacobi equations, for the 
particular case of a particle in a box described by an initially uniform wave- 
function, is connected with the known fractal structure of this wavefunction 
in almost all spatial and temporal directions T]. In section 5 it is conjec- 
tured that, more generally, the nonexistence of these equations corresponds 
to either the wavefunction or its spatial derivative having a fractal struc- 
ture (where the latter case corresponds to examples having finite average 
energies) . 

Conclusions are given in section 6. 

2 A subtlety of the Schrodinger equation 

Attention will primarily be restricted to quantum systems comprising a single 
nonrelativistic spin-zero particle. The corresponding Hilbert space is then 
given by the set of square-integrable complex functions on the configuration 
space of the particle. It is typically assumed, in what follows, that the 
configuration space is one-dimensional (results can easily be generalised to 
higher dimensions). 

Consider first a system that has a Haniiltonian operator H with a dis- 
crete spectrum {-E„}, and corresponding normalised eigenfunctions {ipn^x)} 
satisfying 



A general state of the system at any time t is then specified by 

iP{x, t)=J2 Cne-*^"*/Vn(a;) = e-*^*/V(a;, 0), (2) 

n 

where the coefficients c„ are any set of complex numbers satisfying the nor- 
malisation condition 

E|Cnr = l. (3) 

n 

It follows immediately from Eqs. (Q) and (j21) that one has the identity 

'h ~in{d/dt)]^{x,t) = (4) 

for all states of the system (in particular, one may apply the operator in 
square brackets to each term of the summation to obtain the result). This 
equation is, of course, the Schrodinger equation for the system. However, 
one cannot in general rewrite Eq. (0} in the more familiar form 

in{d/dt)ilj{x, t) = HiPix, t). (5) 

It is this somewhat subtle point, the inequivalence of Eqs. (JH) and (0) for 
certain states, that underlies the main results of this paper. 

As a simple example, consider the case of a particle of mass m confined 
to a one-dimensional box. If the particle is confined to the interval [0,L], 
with H = —h /{2m){d/dx)'^ and the usual (Dirichlet) boundary conditions 
%l){x,t) = at X = and x = L, then the energy eigenfunctions and eigen- 
values are well known to be given by 

ipnix) = (2/L)^/2 gjj^ nTTx/L, En = (nvr/i) V(2mL2), n = 1, 2, 3, . . . . (6) 

For the particular case where the wavefunction is initially uniform over the 
interior of the box, one then has 

dx4'^{x)ip{x,0) = , n= 1,3,5,..., (7) 

with c„ = for n = 2, 4, 6, . . .. Hence, at time t = 0, 

Hi;{x,0) = J2cnEnMx) = —-^J2('^k + l)sm^ -^, 



which diverges for all x G (0,L). Thus, Eq. (0) is meaningless for this 
example: the operator H acts not only to kick the wavefunction out of the 
Hilbert space, but to knock it right out of the set of functions altogether. 

An analogous example of the inequivalence of Eqs. (JD) and ©, for the 
case of a continuous energy spectrum, is provided by a one-dimensional free 
particle of mass m initially confined to some interval, i.e., with Hamiltonian 
operator H = —h'^{d/dxY/{2m) and initial wavefunction 

7/'(a;, 0) = L-^/2gipox/n^ - L/2 < X < L/2 (8) 

(corresponding, for example, to a plane wave incident on a one- dimensional 
slit). It follows that the wavefunction at any later time has the Fourier 
decomposition 

i,{x^t) = (^Y' f dpMP^Po}L/2h^.,./n-^,Hn2nm) 
\TiL) J p-po 

and hence that Htp{x,t) is not well-defined (in particular, the Fourier inte- 
grand of this quantity scales as \p\ for large \p\). 

Thus, for some states, Eqs. (JH) and © are not equivalent - indeed, the 
latter equation has no meaning for these states. It is this fact that lies 
behind the incompleteness of trajectory-based interpretations, as will be seen 
in section 4. First, however, this subtlety of the Schrodinger equation will be 
investigated a little further, in the following section. 

Finally, it is of interest to note that the above examples of inequiva- 
lence arise with respect to the position representation of the quantum state, 
which is of course the representation having fundamental physical significance 
in trajectory-based interpretations. In contrast, no analogous inequivalence 
arises for the Schrodinger equation in the energy and the momentum repre- 
sentations, for either of the examples given above. Moreover, it may be noted 
from Eq. Q that, even in the position representation, the action of the uni- 
tary evolution operator U{t) = e"*-^*/^ is always well defined, even though 
the action of the Hamiltonian H is not. Thus, as per Wigner's theorem jH] 
(and forming a basic element in most axiomatic approaches to quantum me- 
chanics), it is unitary evolution which is fundamental to describing evolution 
on Hilbert space, with the Schrodinger equation following as a secondary 
consequence. 



3 Energy considerations 

It is not difficult to see that the average energy (H) of the above two ex- 
amples is infinite. Thus, while these examples are perfectly valid quantum 
states, it is difficult to conceive of any method for their physical preparation. 
Any interpretation that fails to explain them therefore suffers from a formal 
rather than a physical incompleteness. Hence it is important to explore the 
issue of energy requirements further, and in particular to determine whether 
examples having finite average energies exist. 

It is convenient for this purpose to return to the case of a discrete energy 
spectrum (similar considerations apply to the continuous case), and suppose 
that the amplitudes, energy eigenvalues, and eigenf unctions scale respectively 
as 

|cj~n-", ^n~n^ iV'nl-n'^ (9) 

for large n, with a > 1/2 (to ensure that the state is square-integrable). For 
the examples in section 2 one has a = 1, (3 = 2 and 7 = 0. 

It follows that (ignoring unimportant phase factors in the first line) 



Hip = Yl CnEnlpn ~ Yl 
n n 



n 



/3-a+7 



f3-2a 

n' 



Hence, one can arrange for Hip diverge almost everywhere, while keeping the 
average energy {H) finite, by choosing /3 — a + 7 > and /3 — 2a < — 1 
respectively. This is equivalent to the condition 

(l + /5)/2<a</3 + 7 (10) 

on a (where for consistency one requires that /5 > 1 — 27). Eqs. © and ()10p 
provide a large parameter range corresponding to examples of particles with 
finite average energies for which Eq. (0) is generally invalid. 

Note that the form of the Schrodinger equation in Eq. (0) must necessarily 
be valid whenever Htp happens to be a member of the Hilbert space in 
question (and hence is well defined), i.e., whenever 



(H^) = fdx\HiP\^ <oo. 



Thus all counterexamples must have an infinite expectation value for the 
square of the energy of the system. It follows immediately from Eqs. Q 
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that 13 — a < —1/2, leading via Eq. (fTUI) to the condition 7 < 1/2 for the 
asymptotic scahng of discrete energy eigenf unctions for any counterexample. 
Finally, to give an optical example, consider a single-mode field of fre- 
quency a; in a nonlinear Kerr medium, with photon annihilation operator a, 
number operator A^ = a)a, and Hamiltonian operator 

The energy eigenfunctions ^„(a;), in the usual quadrature representation de- 
fined by X = (a + a^)/2, are Hermite-Gaussians, and for large photon num- 
bers scale as n~^/'^ on any finite interval 9]. Further, for this case one has 
En ^ n^. Thus /? = 2 and 7 = -1/4. It then follows from Eq. ^ that 
amplitudes scaling as 

\cn\ ~ ri"", 3/2 < a< 7/4, (11) 

yield states with finite average energy, for which the form of the Schrodinger 
equation in Eq. (0) is not valid. As a particular example, one may choose 
the initial state of the field to have the number state expansion 

00 
1^0) := [C(13/4)]-V2 ^(n + l)-i3/V), (12) 

n=0 

corresponding to a = 13/8, where C,{z) denotes the Riemann zeta-function 
and \n) denotes the n-th photon number eigenstate. It would be of interest 
to find a scheme for the physical generation of such states. 



4 Incompleteness 

Given any wavefunction ip{x^t) associated with some quantum system, one 
can define quantities P(x, t) and S{x, t) via the polar decomposition 

of the wavefunction. Note that P{x,t) = \ip{x,t)\'^ is the probability density 
associated with finding the system at position x in configuration space, at 
time t. 

In trajectory-based interpretations of quantum mechanics, the position 
of the system is assumed to be a 'real' property at all times, and the prob- 
ability density P{x,t) reflects incomplete knowledge of this property (in the 



de Broglie-Bohm interpretation P{x, t) also has a more fundamental role as 
a 'real' physical degree of freedom associated with the guiding wave 'ip{x, t)). 
Thus quantum mechanics is interpreted as describing an ensemble of systems, 
with each member of the ensemble following a specific trajectory in config- 
uration space. Such a combination of trajectories and statistics provides an 
illuminating quasi-classical picture of quantum systems (although it should 
be noted that the trajectories can have rather 'surrealistic' properties, that 
conflict with naive classical notions of position and its measurement |lUj). 
The aim of such interpretations is to explain the evolution of the ensemble 
as a consequence of the evolution of P and 5". It is necessary to do this for 
all possible wavefunctions if such interpretations are to provide a complete 
explanation of quantum systems. 

Now, if one assumes that the form of the Schrodinger equation in Eq. (jS)) 
is valid, then for the Hamiltonian operator H = —{h /2m)V^ -|- V{x) one 
may multiply this equation on the left by ip* , and take real and imaginary 
parts, to obtain the corresponding equations of motion 

for P and S. The first equation is a continuity equation, ensuring conserva- 
tion of probability, and the second equation is a modified Hamilton- Jacobi 
equation. Both the de Broglie-Bohm interpretation and Nelson's stochastic 
interpretation are based on these equations [H El El EI , and hence the con- 
sistency of these interpretations with quantum mechanics follows whenever 
Eq. (0) is valid. 

However, as was demonstrated by explicit example in the previous sec- 
tions, there are perfectly well defined quantum states for which Eq. (j3)) is 
not valid. For these states one cannot follow the above procedure to de- 
rive corresponding continuity and modified Hamilton- Jacobi equations. It 
follows that interpretations relying on these equations cannot explain the 
evolution of such states, and so are incomplete. This result in fact applies 
not only to trajectory-based interpretations, but to any interpretation based 
on Eqs. p4j) . including the hydrodynamic and exact uncertainty interpreta- 
tions ^i5,,6j. 

Note that there seems at first to be a possible caveat on the above in- 
completeness result. One could argue in particular that such interpretations 
don't need a corresponding Schrodinger equation - they only need Eqs. (fT^ . 
However, such an argument is consistent if and only if Eqs. (J14j) lead to 
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the same predictions as the standard quantum formahsm, for the states in 
question. Unfortunately, this is not the case. In particular, if the continuity 
and modified Hamilton- Jacobi equations are assumed to be a priori valid for 
such states, then the spatial and temporal derivatives of P^/^ and e*'^/^ must 
exist almost everywhere, and one can then derive Eq. (jS)) from Eqs. (fT^ and 
(I14|) . This contradicts the examples of the previous sections. It is therefore 
concluded that the continuity and modified Hamilton- Jacobi equations do 
not correctly describe the states in question. 

For the first example discussed in section 2, of a particle confined to a 
one-dimensional box with the wavefunction initially uniform over the interior 
of the box, the incompleteness of trajectory-based interpretations may be 
seen even more directly. In particular, in such interpretations the initial 
wavefunction corresponds to an ensemble of particles with initial positions 
uniformly spread over the interior of the box (with no particles located at 
the boundaries of the box). Thus, in the neighbourhood of every member 
of the ensemble, P and 5* are initally constant, implying that their spatial 
derivatives vanish. It then follows immediately from Eqs. (jl4|) that P and 
S must remain constant everywhere in the interior of the box - i.e., that the 
ensemble is stationary. This contradicts the quantum evolution, where the 
wavefunction ip evolves as per Eq. (j21). 

The above example further provides an instance of the breakdown of the 
velocity equation 



m ^ 



VS = ih/m) Im {^"V?/'} , (15) 



postulated for each trajectory in the de Broglie-Bohm interpretation P3I2], 
and for the average drift velocity of the forward and backward processes in 
Nelson's stochastic interpretation Ej. In particular, it has been shown 
by Berry that the spacetime dependence of the wavefunction %l){x,t) for this 
example has a fractal structure f7|. Hence Vip is not defined for almost all 
X and t, and therefore Eq. (J15j) cannot be used to define any corresponding 
trajectories or processes. 

5 Fractal connections 

The formal cause of the incompleteness of trajectory-based (and other) in- 
terpretations is seen to stem from the fact that the form of the Schrodinger 
equation in Eq. (0) is not valid for all states - each side of this equation 



can be strongly divergent in the position representation, even for states with 
finite average energy. Here some evidence is collected suggesting that, for 
quantum particles, this divergence is associated with fractal structures of the 
corresponding wavef unctions. 

For the example of the one-dimensional particle in a box, with a wave- 
function initially uniform over the box. Berry has shown that the probability 
distribution P{x,t) has fractal dimension 3/2 in the spatial direction for al- 
most all fixed times t, and fractal dimension 7/4 in the time direction for 
almost all fixed positions x. A simplified expression for P{x,t) for this ex- 
ample is given in Ref. ^T], and an approximate experimental realisation of 
the fractal structure, via an optical analogue, is discussed in Refs. jT2l CH] • 
Further fractal and near-fractal wavefunctions have been constructed by Wo- 
jcik et al. (U] and by Amanatidis et al. [T^ (all having infinite average energies 
in the fractal limit). For such wavefunctions both P and S in Eq. (fT^ are 
typically also fractals, and hence provide further examples where Eqs. fjl4|] 
and (fT3j) are not well defined . 

The fractal nature of the wavefunctions in Refs. [3 E] was derived as a 
consequence of the result that functions of the form 

/(x) = ^a„e^«^ 

n 

for which the amplitudes a^ scale asympotically as 

|a„| ~ |n|"'= with 1/2 < z < 3/2, 

are continuous but nondifferentiable, and have fractal dimension 5/2 — ^ [7]. 
Thus, for example, for a particle in a one- dimensional box with amplitudes 
\cn\ ~ 'n~"' in Eq. Q, it follows via Eq. (jU)) that the wavefunction has a 
fractal structure whenever 1/2 < a < 3/2 (this includes the particular case 
of the initially uniform wavefunction, for which a = 1). Note, however, that 
since (H) ~ Z]n^^~^°; ^-U such examples have infinite average energies and 
hence cannot be physically prepared. 

It turns out that fractal structures can also be associated with states hav- 
ing finite average energies. In particular, consider again states of the particle 
in the box with coefficients \cn\ ~ n"" in Eq. (0), where a is now chosen 
to be in the parameter range corresponding to Eq. (fTUj) in section 3. Since 
(3 = 2 and 7 = for this case, this range is given by 3/2 < a <2. Now, by 
construction, the corresponding states have finite average energy and do not 
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satisfy the form of the Schrodinger equation in Eq. (0). On the other hand, 
they do not satisfy the fractal criterion given above. However, from Eq. (jUj) 
one finds that the spatial derivative of the corresponding wavefunctions has 
the form 

{d/dx)il){x,t) = v2 TiL^^'"^ y^ ncn cos m:x/ L. 

n 

This quantity does satisfy the above fractal criterion whenever 1/2 < a — 1 < 
3/2, and hence in particular for the range 3/2 < a < 2 of interest. The 
corresponding fractal dimension is 7/2 — a, and hence also lies between 3/2 
and 2. 

Based on the above results, it is conjectured that the incompleteness 
of trajectory-based interpretations for quantum particles corresponds to the 
existence of states for which the spacetime dependence of the wavefunction, 
or of its spatial derivative, has a fractal structure. 

6 Conclusions 

The incompleteness of trajectory-based interpretations arises for systems 
with unbounded Hamiltonian operators. It applies not only to the very sim- 
ple case of a wavefunction initially uniform over the interior of a box, but also 
to a number of examples having finite average energies and hence which can, 
in principle, be physically prepared. For quantum particles the incomplete- 
ness of such interpretations appears to be connected with associated fractal 
structures. 

Strictly speaking, one should differentiate here between the notions of 
formal and physical incompleteness. The results of the paper show that 
trajectory-based interpretations are formally incomplete, as they do not de- 
scribe all possible states in the Hilbert space, even though these states and 
their (unitary) evolution are well defined. The results suggest that such in- 
terpretations are also physically incomplete, as some of the states in question 
have finite average energies, and so can plausibly be physically prepared from 
finite resources. However, it is open to proponents of such interpretations to 
argue for physical completeness on the grounds that even the finite-average- 
energy counterexamples are unphysical. For example, noting the discussion 
in section 3, it would suffice to provide a convincing argument that all mo- 
ments of the energy of a physical system must be finite. 

In the case of wavefunctions for which P and S are fractals there would 
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appear to be little chance of overcoming incompleteness, via some supple- 
mentary rule that specifies how to generate the corresponding (presumably 
fractal) trajectories. In contrast, in the case of wavef unctions for which the 
spatial derivative is a fractal, the velocity equation in Eq. (fT^ is well de- 
fined, and hence might be used to specify a set of associated trajectories 
(eg, via i = f in the de Broglie-Bohm interpretation). However, given the 
nonexistence of a corresponding continuity equation as per Eq. (I14|) . it is not 
clear that an ensemble of these trajectories can evolve in agreement with the 
Schrodinger equation in Eq. (JD). It would be of interest to perform some 
numerical experiments in this regard. 

As previously remarked, the incompleteness result applies to any interpre- 
tation that relies on the continuity and modified Hamilton- Jacobi equations 
in Eq. (fT^ [H I2l El IH ISl E]- Thus quantum mechanics goes where these 
intepretations do not follow, despite their (at least in principle) duty to do 
so. 
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